Spin distribution of nuclear levels using static path approximation with random-phase 

approximation 
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We present a thermal and quantum-mechanical treatment of nuclear rotation using the formalism 
of static path approximation (SPA) plus random-phase approximation (RPA). Naive perturbation 
theory fails because of the presence of zero-frequency modes due to dynamical symmetry breaking. 
Such modes lead to infrared divergences. We show that composite zero-frequency excitations are 
properly treated within the collective coordinate method. The resulting perturbation theory is free 
from infrared divergences. Without the assumption of individual random spin vectors, we derive 
microscopically the spin distribution of the level density. The moment of inertia is thereby related 
to the spin-cutoff parameter in the usual way. Explicit calculations are performed for 56 Fe; various 
thermal properties are discussed. In particular, we demonstrate that the increase of the moment of 
inertia with increasing temperature is correlated with the suppression of pairing correlations. 
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I. INTRODUCTION 

The spin distribution of nuclear levels is important for 
Hauser-Feshbach-type calculations of astrophysical reac- 
tion rates [![. Generally, the spin distribution of nuclear 
levels is assumed to be given by 



p{EJ) = p(E) 
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where p(E, I) is the level density for a given spin I and 
a is the spin-cutoff parameter 0, H, 0, H, • It should 
be noted that the total level density (summed over all 
spins) is given by p(E) = J2j p(E, I), while the total 
state density given by W{E) — £j(2/ + l)p(E,I) con- 
tains an additional factor 21+1 for each level in order to 
take into account the m degeneracy. In deriving the spin 
distribution, one typically assumes that the individual 
nucleon spins are pointing in random directions, hence 
the spin distribution becomes Gaussian and can be de- 
scribed by only one parameter, the spin cutoff parameter. 
This parameter can be related to the moment of inertia 
frigid by 



(3) 



where T is the thermodynamical temperature This 
result can be derived within an independent-particle 
model with individual nucleon spins but neglecting cor- 
relations Q. Microscopic calculations of the spin distri- 
bution are difficult when correlations are present. Exper- 
imental data on the spin-cutoff parameter are available 
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only in a few cases [4|. Therefore, theoretical predictions 
of the spin-cutoff parameter are greatly desired. Alhas- 
sid et al. have recently investigated the spin distribu- 
tion in the framework of the shell model Monte Carlo 
(SMMC) M, Q and the static path approximation (SPA) 
method .8). Assuming that the spin distribution can be 
described by Eq. they found a significant suppression 
of the moment of inertia and an odd-even staggering of 
the spin-cutoff parameter due to pairing correlations at 
low excitation energies. 

The SPA method @, [13, d El is a useful treatment 
to evaluate approximately the partition function of finite 
systems with separable interactions such as the pairing 
plus quadrupole-quadrupole (P+QQ) interaction, and it 
provides the exact result at high temperatures. However, 
as temperatures decrease the SPA becomes inaccurate 
because quantal fluctuations around the mean field can 
no longer be neglected. In leading order perturbation 
theory, small amplitude fluctuations give corrections to 
the partition function which can be treated by random- 
phase approximation (RPA) [H [3 GJ, Qj] . 

However, such a perturbative treatment exhibits an in- 
frared problem in the case where the nucleus is modeled 
in terms of independent nucleons moving in a deformed 
mean field: there are zero-frequency modes among nor- 
mal vibrational modes. This can be understood because 
the SPA method breaks the symmetry of rotational in- 
variance for a deformed nucleus. It operates in a body- 
fixed reference frame which undergoes time-dependent 
transformations. The collective (zero-frequency) excita- 
tion, i.e., the motion of the rotating frame, is treated on 
an equal footing with other excitations. This leads to a 
Lagrangian which displays local symmetry yielding the 
mechanical equivalent of gauge invariance. The collec- 
tive coordinate method developed by Gervais and Sakita 
[l7| treats the inherent symmetries of the problem consis- 
tently. It separates the so-called spurious motions from 
the intrinsic vibrations. In the laboratory frame, these 
zero-frequency excitations are not really spurious; they 
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represent a rotational motion which has to be treated 
separately. 

In the present paper, we describe microscopically the 
nuclear rotation of a deformed nucleus including thermal 
and quantal fluctuations using the SPA+RPA method, 
and derive the associated spin distribution of nuclear lev- 
els. First, using the collective coordinate method, the ro- 
tation is described in terms of the zero-frequency modes 
which emerge from the breaking of rotational invariance 
of the deformed mean field. Secondly, the spin distri- 
bution of nuclear levels is calculated by a proper treat- 
ment of the zero-frequency modes in the framework of 
the SPA+RPA method. In Sect. HH the collective coor- 
dinate method is applied to the partition function of an 
interacting fermion system. In Sect. IIHI we show how 
to calculate the partition function of the P+QQ model 
using the SPA+RPA method. In Sect. IIV1 we present a 
correct treatment of the zero- frequency modes and derive 
the spin-dependent level density. In Sect. [V] the model 
is applied to 56 Fe and related thermal properties are in- 
vestigated. Concluding remarks are given in Sect. I VII 



II. THE COLLECTIVE COORDINATE 
METHOD 



According to the theory developed by Gervais, Jevicki, 
and Sakita [l8j |. a problem involving only fermion de- 
grees of freedom can be transformed into an equivalent, 
constraint problem including both fermion and collective 
degrees of freedom. Concerning the problem at hand, 
a perturbative treatment of nuclear rotation within the 
path integral formalism has been proposed using the col- 
lective coordinate method [HI, [2(|. Let us consider a 
partition function of an interacting fermion system given 
by the path integral after Wick rotation 



Z = I ]JV[cl]V[c a ]t 

J a 

Here, the action Se{P) is given by 



-S B (/3) 



(4) 




■H{c\c) dr, 



(5) 



where (3 = 1/T is the inverse temperature and the dot 
indicates a derivative with respect to the imaginary time 
r. The antiperiodic boundary condition c{(3) = — c(0) 
arises since the operators c are Grassmann numbers due 
to the fermionic nature of the field. The action Se{0) is 
rotationally invariant with respect to the SO (3) algebra. 

We shall now be concerned with the case in which the 
fermion system becomes deformed and breaks the SO(3) 
invariance. As mentioned in the Introduction, a naive 
perturbative treatment fails due to the presence of zero- 
frequency modes. A better way to tackle this problem is 
to introduce collective coordinates. We therefore adopt 
the path integral formulation of Gervais and Sakita [17| . 



In order to separate the collective coordinates, we employ 
the equality 

i= fWviuW^hU-j^.b)) 

x,5(a ! (fe t ,6))det[J,e]pB, (6) 

where 0^ (b ' , b) is an arbitrary function of the fermion 
fields (b' , b) with the condition that the determinant of 
the Poisson bracket [ ]pb does not vanish. The fermion 
fields in the body- fixed frame (w , b) are obtained by a 
linear transformation from the laboratory frame 



i ) = p R (<i>i)c a p R (ct> i )- 1 , 



(7) 



where -Pa(<^i) is the rotational operator with the Euler 
angles (pi (i = 1,2,3), and Li are the three components 
of the collective angular momentum. 

Inserting the Eq. © into the right hand side of Eq. ((31) 
and exploiting the constraint fiLi = Ji (b^ , b) , one finds 

Z = [ l[V[bt]V[b a ] I Y[V[Li]V[<lH] 

x exp f-5 B (/3) + ^ ^n j (r)^(T)drj 

xS{hL z ~ Ji(b\ b))5(Qi(tf, b)) det(y) det[J, 0] PB) (8) 

where IT are conjugates of the angular variables 4>i de- 
fined by the linear transformation 



— sin 4>2 cos (f>3 sin cf>3 
I' = | sin (j>2 sin 3 cos fa 
cos 4>2 1 



(9) 
(10) 



Let us now define Zj by imposing on Eq. © the bound- 
ary condition Li — * Ii for r — ► or j3, where the total 
angular momentum I is given by I — Then 
the path integral over 4>i and Li can be computed by 
standard methods: 

Zi = J d*I J Y[v[bi]V[b a ]e- s ^5 (i- ] JY^\ 

xS(hh - Ji(tf,b))6(Qi(bt,b)) det[J, 8]p B . (11) 

It is now useful to specify the <di(b\b) as the conjugate 
angular variables of the Ji{b\ b). Then, the gauge condi- 
tion 8(Qi(b\b)) fixes the position of the intrinsic frame 
of reference relative to the rotating body. This eliminates 
rotations of the system as degrees of freedom associated 
with the fermion fields. However, the components Ii of 
angular momentum in the intrinsic frame depend upon 
the choice of the gauge angular variables <di(b\ b). 
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III. MODEL 

A. The P+QQ Hamiltonian 

Let us consider the P+QQ Hamiltonian [2l|, [22| as a 
model Hamiltonian 

H = ff s . p . - GPtp _f£ Q^Oap, (12) 



a 
a 

Q 2M = ^(a|r 2 y 2/ ,|^4^, 



(13) 
(14) 
(15) 



where i? s . p . is a single-particle Hamiltonian, P is the 
monopole pairing operator, G is the strength of the 
pairing interaction, are components of the mass 

quadrupole tensor, and x is t ne strength of the quadru- 
pole-quadrupole interaction. This model has been ap- 
plied to the SPA+RPA method including thermal and 
quantal fluctuations [23| . 



B. The Hubbard-Stratonovich transformation 

Using the Hubbard-Stratonovich transformation [24 
const = f VCVCY[V[a;}V[a^} 

[\tG(C - P + )(C - P) 
Jo 



x exp 



x exp 



2/W 



,(16) 



the partition function of Eq. (1111) can be written as the 
auxiliary field path integral [la, [25| 



Z I = / d 3 I / VCVC / n^KPM T r 



-fJH' 



x exp 



, , Kl 2 x, ,a 



xS{Mi - Ji(b\ b))5{Qi{b\ b)) det[J, 9] PB , 
where 



Tr 



-/3H' 



= /n^aMMexp 



(17) 



(18) 



C. The SPA+RPA method 

In the perturbative expansion of the partition function 
Zj around the static fields, the zeroth-order term gives 
the Hartree-like solution, while the second-order terms 
lead to the RPA corrections. This approximation may be 
largely improved by applying it to all time-independent 
paths of C an d c M . We therefore expand the partition 
function (|17p around the static paths £ and in Fourier 
series, 



(20) 
(21) 



n#0 



where u n = 2nn/ (3 are the Matsubara frequencies. Let 
us now define the principal axis of the quadrupole poten- 
tial in the intrinsic frame by <7o = fi.wo/?2 COS7, <t 2 = 
(T_2 = fuo§{3i sin7/v / 2, and o\ — er_i = 0, where 
Hu)q — 41 MeV /A 1 / 3 . The pairing fields can be rewritten 
as £ = Ae - ^ using real values A and i/j. 

It is now convenient to introduce quasiparticles by di- 
agonalizing the Hamiltonian of Eq. (|19p 



w v-t 



(22) 



where the matrix W satisfies the unitarity condition 
yytyy = 1. According to the Bloch-Messiah theorem 
[1^], the above unitary matrix can be decomposed into 
three matrices 



w = 





) 








1 





D* j 




I) 





c* ) 



(23) 



The first transformation D is determined by diagonaliz- 
ing the deformed term of Eq. (TIT)]) 

Hdef = Hs.p. ~ X^ ^nQiii 



H, 



s.p. 



f^l I Q22 + Q2-2 . \ , .N 
x ^Q 2 o cos 7 H -= sin 7J , (24) 

from which we obtain the deformed fcrmion fields 

d k = J2 D k,*b a , (25) 

a 

where D ka are the matrix elements obtained by solv- 
ing the eigenvalue problem H^\k) = e k \k). Now, the 
Hamiltonian H' of Eq. (|19[) can be simplified to 



H' = J2^k4d k -A(P^ + P) 



(26) 



The matrices U and U of Eq. (f2"3"| diagonalize the pairing 
term of Eq. (|26|) . The matrices are diagonal and deter- 
mined by solving the Hartree-Fock-Bogolyubov equations 



£ k A 
A -e k 



Uk 
Vk 



= E k 



Uk 

v k 



(27) 
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where E k are the quasiparticle energies. 

In a small system such as the nucleus in which the 
particle number is strictly fixed, the canonical partition 
function should be employed. However, exact number 
projection is difficult to implement in the SPA+RPA 
formalism; it is easier to apply number parity projec- 
tion. Performing the Gaussian integral over and T]^ n 
in second order perturbation theory and introducing the 
number parity projection P s = (l + se' l7vN )/2 where s de- 
notes the even or odd number parity (l5| , we obtain the 
partition function 



G \2n J 

poo 

dAA I dfafe 



/O JO ./n 

'A 2 



r/3 



dy \ sin(37)| 



x 4ttI exp 



1 



G 



+ 0^2 P 



Tr 



P s e 



-0H' 



xS(hI t - J 4 ( a t, a )) ( 5(a t (at, a ))det[J,e]pBC R PA, (28) 

where the parameter k = (huo/b 2 ) 2 j\ and b oc A 1 / 3 
is the harmonic oscillator length. The expression l|28p 
contains both the static mean-field contributions and the 
associated Gaussian corrections. These small-amplitude 
quantal corrections are included in the factor 



Crpa = Crpa(Q)Crpa(A), 



(29) 



where Crpa(Q) and Crpa (A) are the corrections due to 
the QQ and pairing correlations, respectively. The SPA 
partition function is obtained by neglecting these RPA 
corrections. Crpa(C) is explicitly given by 



Crpa(Q) = n det (! ~ x-RK)r x . 



(30) 



n>0 



Here R(to n ) is the response function matrix given by 
QkiQkiiEk+EOil-fk-h) 



i?K) = ]T 



[Ek+E^+ul 

)klQkl{Ek — Ei)(fk — fl) 



Qm = qki(ukvi + v k ui), 
Qu = qu{ukui - vkvi), 



(31) 

(32) 
(33) 



where q k i are quadrupole matrix elements and f k are 
Fermi occupation probabilities f k = (1 + e l3Ek )^ 1 . The 
response function can be calculated by solving the disper- 
sion equation det(l — xRi^n)) = 0. It can be also ob- 
tained by diagonalizing the finite-temperature RPA equa- 
tions 0, mi 

{[H WK (Q),Ol]) = m v {Ol), (34) 
where f2„ are the RPA frequencies and the RPA modes 



are defined by 

kl 

J2( Z M a k ai ~ Zkl a W)- 



(35) 



kl 



The symbol ( ) in Eq. |34|) denotes thermal averaging 
with respect to quasiparticles according to Wick's theo- 
rem 



(F) = Tr 



Tr 



Fexp ^-(3 E k a\.ak 
exp l-P^E* 



Z qp , (36) 



(37) 



where, as an an example, we give (atai) = S k if k - At 
this point, it would be more correct to use thermal aver- 
age and Fermi occupation probabilities in number parity 
projection instead, but the corresponding corrections in 
the RPA are normally small, and would not change our 
conclusions. It would be important for the odd-even ef- 
fects in the moment of inertia [|[ . The RPA Hamiltonian 
^rpa(Q) can now be expressed in a diagonal form 



(38) 



There are, however, zero-frequency modes O„ = o due to 
the symmetry breaking of the deformed Hartree solutions 
on the static path. The correct treatment of these modes 
will be given in the next section. The RPA correction 
Crpa(Q) can now be written exactly in terms of the RPA 
frequencies VL V 



Crpa(Q) 



kl 



sinh 



2 



iLlfsinh^ 



(39) 



where the prime in Y[' k i denotes the restriction of the 
product to pairs (k, I) that satisfy the conditions k < I 
and E k + Ei 7^ 0. Note that for deformed or heavy nuclei 
there are many RPA frequencies Vt v in the numerator and 
many E k + Ei pairs in the denominator. 

The pairing RPA correction Crpa (A) is obtained in a 
similar way as Crpa(Q) by 



Crpa (A) = JJ 



ujk sinh[f3E k ] 
2Eksmh[/3cok/2)' 



(40) 



Here, uj k are the conventional thermal RPA energies and 
E k = + A 2 with e' k = Ek -n- G/2. 



IV. ROTATIONAL MOTION AND 
ZERO-FREQUENCY MODES 

We consider the two-body Hamiltonian H which is in- 
variant under a continuous rotational operation gener- 
ated by the angular momentum operators Ji. However, 
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the deformed solutions on the static path violate this 
symmetry and lead to zero-frequency modes which are 
regarded as spurious modes in the intrinsic frame. Then, 
the rotational invariance must be restored by the residual 
interaction, which is defined as the difference between the 
exact and the mean-field Hamiltonian. Proper inclusion 
of the residual interaction will therefore restore the rota- 
tional invariance and will provide the rotational energy. 
To achieve this goal, we present in this section the cor- 
rect treatment of zero-frequency modes within the finite 
temperature RPA. It is shown that the zero-frequency 
modes in the intrinsic frame are not spurious in the lab- 
oratory frame and correspond to three-dimensional rota- 
tional motions. 



A. Treatment of zero-frequency modes 

For the calculation of the partition function of the con- 
straint system (see Eq. (|28p). we follow the treatment 
proposed by Marshalek and Weneser [29| ■ For the sake of 
simplicity, we assume an axially symmetric mean field for 
the deformed Hartree calculation with one violated sym- 
metry and neglect the 7 degree of freedom, for simplicity. 
In particular, we consider the case where the symmetry 
is violated for rotations around the x axis which is per- 
pendicular to the symmetry or z axis. Three-dimensional 
rotation can be treated in a similar way. 

When the eigenvalue equation has zero-frequency 
solutions due to the breaking of rotational invariance in 
the intrinsic frame, 1 the usual treatment is known to 
cause problems concerning the completeness and normal- 
ization. Therefore, instead of O v and Oj, we define the 
following coordinate and conjugate momenta 



Vv 



^(O v +Ot), (41) 

(42) 



The finite-temperature RPA equations are then ex- 
pressed as 



([Hrpa, <Zv]pb) 
([Hrpa, Pv]pb) 



= ihSll(q u ). 



(43) 
(44) 



For the zero-frequency mode O„ = o in the finite-tempera- 
ture RPA equations (|34|) . the angular momentum J x and 
the canonical conjugate angle coordinate Q x are defined 



by J x 



Jx /2 po and Q x 



Jx 1 qo, respectively, where 



J x is the moment of inertia. and J x can be deter- 



1 Note that the zero-frequency mode (Goldstone mode) in the 
mean field plus RPA may become non-zero (imaginary or com- 
plex) in the SPA+RPA. 



mined by solving the equations of motion 

([ITBPA.ejpB) = ^(Jx), (45) 

{[Qx,Jx]pb) = ih, (46) 

where Eq. f?5] simply corresponds to the self-consistency 
condition Q x = [#rpa - J%/2Jx, ©x]pb = 0, and we 
impose the gauge fixing condition <d x =0 of Eq. f28|) . 
-Hrpa(Q) can now be separated into an intrinsic and ro- 
tational part 

ffRPA(O) = J2 m „OtO u + ^r. (47) 

Thus, the zero-frequency mode can be eliminated from 
the RPA equation and Crpa(Q) can be written as 



Crpa(Q) = C^ PA (Q)e-* 



(48) 



where Cp PA (Q) is the correlation factor in the absence 
of zero-frequency modes. 

In a similar fashion, spurious modes should be elimi- 
nated from the pairing RPA corrections of Eq. (|4"0"1) . This 
is not done in the present work since (i) exact number 
projection is needed for this procedure and (ii) the cor- 
rections are thought to be small in our number parity 
projection scheme. 



B. Moment of inertia 

In order to connect with the result of Thouless and 
Valatin [3(|, we introduce the function G 



G = iJ x Qx/n = Yl 
k,l 



UkVia\a\ - Vkuiajak) 



- 9u{ukUia\ai - v k via\ak) ■ (49) 
Then the Thouless- Valatin equations become 

([H RP a,G]pb) = (Jx), (50) 



([Jx,G]pb) — Jx 



(51) 



Inserting Eq. (|49|) in the Thouless- Valatin equations 
yields 



ft (jx) k l = [ E k + Ei)9kl - [QkWmn ~ Q 



-Qkig*mn ~ QklQ 



mil 



kl9n 



Qrnn, (52) 



ft Ux)kl = (Ek — Ei)gkl — X^^l \_Q mn 9mn — Qmn9mn 

ki, (53) 



~\~Qmn9 m n Qmn9 

and the moment of inertia 

Jx = h^2(j x )ki [gikjukVi - Vkujfjl - f k - fi) 



k.l 



+9ik(u k ui + v k vi) 2 {f k - //)] , (54) 



where we used the relations 



9lk 



gik 



~9kl> 
-Vz. 



(55) 



Neglecting the two-body potential matrix elements in 
Eqs. (|52")) and (JHSJ), gik and gi k can be approximated by 



9kl 



(jx)k 



'{Ek + EiY 



9kl 



u 



x )kl 



(E k - Ei) 



(56) 



and the moment of inertia at finite temperature simplifies 
to 



J x = 2H 2 



(E k + Ei) 



(u k vi - v k ui) 2 (l - f k - fl) 



k,l>0 

+2h 2 J2 ^ Jx)k l s (u kUl + v k v{f{f k - - ./, ). (57! 



k,l>0 



(E k — Ei ) 



In the T — ► limit, the moment of inertia becomes the 
well-known Belyaev formula [3l| 



J X = 2K 2 (jx)2hl 



k,l>0 



(Ek + Ei) 



(58) 



(u k vi - V k Uif 



Inserting Eqs. P^ - pS")) into the partition function 
under consideration of the gauge fixing conditions HI = 
J x and Q x — 0, we obtain the partition function for an 
axially symmetric nucleus 

Z s j = |? f^Y ' dAA J™ d(3 2 ^Z s (A,(3 2 ) 



x AttI exp — 



2r2 



h z i 

2J X 



(59) 



where 



Zs(A,p 2 ) 



exp 



'A 2 1 
-G + 2' 

x^ qp ^ PA (Q)^RPA(A), (60) 



+ ^01 ) 



and Zf p is the quasiparticlc partition function with the 
number parity projection 



Z ? = III* 



(1 



l + s]Jtanh 2 (^/3) 



(61) 



with 7 fe = e k - fJ. - E k . 

C. Level density 
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FIG. 1: (Color online) Experimental (squares) and calculated 
level densities p(E) as function of excitation energy in 56 Fe. 
Full circles and triangles denote the level and state density, 
respectively, from SPA+RPA calculations. The solid line is a 
back-shifted Fermi-gas model. 



In the saddle-point approximation, the level density is 
given by 



P(E,I) 



Z s ,ie 



\2-k3 2 InZsj/dp 2 } 1 / 2 ' 



(63) 



We will now apply the saddle-point approximation 2 to 
the two integrals in the partition function Z s j of Eq. 
(|59|) . However, Eq. ((59)) still contains the 21+1 degen- 
eracy connected to the magnetic quantum number. In 
order to determine the level density (as opposed to the 
state density), we make the transformation 7 — > 7 + 1/2 
and J 2 — > /(J+l) in Eq. ([59)1 from a classical to a quantal 
angular momentum and divide by 27 + l. 3 Entering the 
thusly modified partition function (|59| into the saddle- 
point approximation for the level density (f6"3")> . we finally 



The level density for a system with a fixed number of 
particles can be evaluated by an inverse Laplace trans- 
formation of the partition function Z s j 

p(E,I) = — d^ E Z sJ (P). (62) 

Lm JPa-ioo 



2 The maximum in the saddle-point approximation leads to an 
effective mean-field equation, hence, a sharp phase transition as 
predicted by the ordinary mean-field equation is avoided. 

3 In order to keep in contact with experiment, we will consider 
this pseudo partition function for the remainder of this work and 
determine consistently all quantities from there. 
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FIG. 2: (Color online) Heat capacities as a function of tem- 
perature T in 56 Fe. The solid and dotted curves denote the 
calculated heat capacities in the SPA+RPA and SPA, respec- 
tively. Both curves exhibit clear S shapes. 



obtain 



P(E,I) 



1) exp 



h 2 i(i + i) 



where the total level density is given by 



P{E) 



Z s e? E 



[2nd 2 InZs/dp 2 } 1 / 2 



(3 , (64) 



(65) 



with 



5/2 



APIZ S (A,(3 2 ) 



and B is the fluctuation matrix 



Bj 



d 2 ln(A^Z s (A,(3 2 )) 
dyidyj 



(66) 



(67) 



withy = (A,/3 2 ). 

The level density ([64| is the relevant spin-dependent 
microscopic level density for reaction cross-section calcu- 
lations in the presence of correlations. By entering the 
partition function (|60)l in the Eqs. ([66]) and (|67|). we can 
calculate this level density as a function of excitation en- 
ergy, where the thermal energy can be calculated from 
E = — din Z s /dj3. Comparing this level density with the 




T (MeV) 



FIG. 3: (Color online) Effective pairing gaps (a) and their 
derivatives (b) as a function of temperature for 56 Fe. The 
solid and dotted curves are for neutrons and protons, respec- 
tively. 



spin-cutoff model given in Eq. (p}, the spin-cutoff pa- 
rameter a 2 can be related to the moment of inertia J x 
by 



a 2 = 



(68) 



In the next section, we shall perform numerical calcula- 
tions and show that there is a notable difference between 
the shape of the level density and the state density which 
contains the factor 21+1 due to the m-degeneracy. 



V. NUMERICAL CALCULATIONS 

We consider 56 Fe [1,1131 as an illustrative example. Re- 
cently, the level density in 56 Fe has been measured nearly 
up to the neutron binding energy [33l [34], |35| . In our cal- 
culation, we use the single-particle energies £k given by an 
axially deformed Woods-Saxon potential with spin-orbit 
interaction [36J. The Woods-Saxon parameters are cho- 
sen such as to approximately reproduce the experimental 
single-particle energies extracted from the energy levels of 
the odd nucleus 41 Ca (a 40 Ca core plus one neutron). 15 
doubly-degenerate single-particle levels for neutrons and 
protons outside the 40 Ca core are considered; continuum 
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FIG. 4: (Color online) Moment of inertia (a) and the deriva- 
tive (b) as function of temperature for 56 Fe. The dashed line 
shows the rigid-body moment of inertia, the solid curve is the 
result of the SPA+RPA calculation, and the dotted curve is 
a fit to the calculation according to Eq. (|73[) . 



levels are neglected since their contributions are small Q . 
We adjust the pairing- force strengths at G n = 25 /A MeV 
and G p = 29/A MeV for neutrons and protons, respec- 
tively, and the QQ-force strength at x = 240 J A4 5/ ' 3 & 4 . 
As mentioned above, the SPA+RPA breaks down at low 
temperature. However, it has recently been shown that in 
the monopole pairing case the SPA+RPA with number- 
parity projection reproduces well exact results even for 
low temperatures |15| . Therefore, number-parity projec- 
tion is essential to describe the thermal properties at low 
temperature. Calculations are shown in Fig. [TJ The level 
density calculated with the SPA+RPA reproduces fairly 
well the slope and magnitude of the experimental data. 
On the other hand, beside the expected difference in mag- 
nitude, the calculated state density which contains the 
m-degeneracy shows also a different slope on the loga- 
rithmic plot. The calculations thus indicate an increase 
of a oc (I) with excitation energy, see also Eq. ([68]) . 

FigureJT] also shows a back-shifted Fermi-gas model 
SroSH according to 
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FIG. 5: Correlation of the average effective pairing gap with 
the moment of inertia. The straight line is to guide the eye. 



where the back-shifted energy is U = E — E\ and the 
spin-cutoff parameter a is er 2 = 0. 0888 A 2 ' ' 3 yfaU [3]. The 
level-density parameter a and the parameter E\ are given 
by a = A/10 = 5.6 MeV" 1 and E 1 = 1.0 MeV, respec- 
tively. The level density from SPA+RPA calculations 
is close to the back-shifted Fermi-gas model for a wide 
range of excitation energies. 

Figure [5] shows the calculated heat capacity for 56 Fe. 
The SPA+RPA deviates from the SPA result for tem- 
peratures above T =0.9 MeV; quantal fluctuations be- 
come important in this region. This behavior is consis- 
tent with the one found in Ref. [37[ ■ Both of our results 
exhibit a characteristic 5-shape around T c ~ 0.9 MeV. 
This critical temperature is higher than the SMMC re- 
sult 4 ~ 0.7 MeV [II] but it is still smaller than the ex- 
perimental one ~ 1.3 MeV [33| . Level densities in sev- 
eral nuclei have recently been measured [H, H(| and the 
suppression of pairing correlations was studied in detail 
[H S3]. S-shaped heat capacities with T c ~ 0.5 MeV 
were also observed in 162 Dy, 166 Er, and 172 Yb 0,111; 
they were interpreted as a signature of the breaking of 
nucleon Cooper pairs. We suggested that the suppression 
of pairing correlations at finite temperatures should also 
appear in the thermal odd-even mass differences [4l|, \^ . 
In general, the 5-shape of the heat capacity has been 
attributed to the reduction of the pairing energy which 
can be calculated from G T (P f P) = G T Td In Z s /dG T [H. 
In fact, this pairing reduction can be seen in thermal 
odd-even mass differences extracted from the experimen- 
tal level densities [HI, H2, |H|. On the other hand, it 



Pbsfg(E) 



exp 2 V aU 



12V2V/ 4 [/ 5 /V 



(69) 



4 We would like to emphasize that in contrast to the present cal- 
culation and the experiment, the SMMC result is obtained from 
a partition function which includes the m degeneracy 
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has been argued that the 5-shape might be accounted 
for as an effect of the particle-number conservation on 
the quasiparticle excitations j4f|. However, even with- 
out exact particle-number projection, the SPA gives an 
5-shaped heat capacity. In this sense, we emphasize that 
this S shape cannot be explained solely as an effect of 
the particle-number projection. 

Figure^ shows the effective pairing gap defined by 



1 d In Z s 

dG T 



1/2 



(70) 



A™ ff and A^ ff decrease in a similar fashion around T — 
0.9 MeV, although the proton pairing is a little bit 
stronger than the neutron pairing. We can now iden- 
tify the inflection point of the A e fj curves as the critical 
temperature of a pairing phase transition. In order to 
determine the position of the inflection point precisely, 
we differentiate the A e g curves with respect to the tem- 
perature, see Fig.[3jD. From the peaks of these derivatives 
we read off a critical temperature around T = 0.9 MeV. 
Hence, the suppression of A e g is well correlated with the 
S shape of the heat capacity. As pointed out in our previ- 
ous paper (4lj , the critical temperature T c is proportional 
to the pairing gap A. This is expressed by the relation 
T c « 0.57A from Bardeen-Cooper-Schrieffer (BCS) the- 
ory. The critical temperature T c sa 1.0 MeV estimated 
from this relation is close to the position ~ 0.9 MeV of 
the peaks in Fig. [3b. 

We will now apply a similar analysis to the moment of 
inertia. Figure [4^i shows the moment of inertia as func- 
tion of temperature T. For comparison we also show 
the rigid-body value J^ ld = 0.0137A 5 / 3 MeV" 1 de- 
noted by a dashed line. The calculated moment of in- 
ertia is smaller than the rigid-body value at low tem- 
perature. This strong suppression is well known as the 
effect due to the presence of pairing correlations in low- 
lying, superfluid-likc BCS states. For high temperatures, 
however, the calculation approaches the rigid-body value 
because the nucleon Cooper pairs break with increasing 
temperature. Again, we can define the critical tempera- 
ture of the pairing phase transition as the inflection point 
which is determined precisely by differentiation. Differ- 
entiation of the J x curves reveals peaks around the criti- 
cal temperature of T = 0.9 MeV (see Fig. 2b). Hence, the 
increase of J x is well correlated with both the S shape 
of the heat capacity in Fig. [2] and the suppression of the 
effective pairing gaps in Fig. [3J 

To investigate more closely the relation between the 
effective gaps and the moment of inertia, we exam- 
ine the correlation of the average effective pairing gap 
A e g = (A" ff + A^ ff )/2 with the moment of inertia, see 
Fig. [5] Above J x w 5.0 fi. 2 /MeV, which corresponds 
to T i=s 0.7 MeV, this correlation is close to a straight 
line, while deviations from the line become apparent be- 
low this temperature. The observation of such a corre- 
lation across the critical temperature of T w 0.9 MeV 
means that the increase of the moment of inertia can be 
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FIG. 6: (Color online) Spin distribution of level density at 
different temperatures T for 56 Fe. The SPA+RPA results 
(solid curves) are compared to the spin-cutoff model using the 
rigid-body moment of inertia (dashed curves) and the results 
from the fit to the moment of inertia according to Eq. f|T3p 
(dotted curves). The fit and the SPA+RPA calculations are 
almost indistinguishable. 



attributed to the suppression of the pairing correlation. 
This conclusion has already been obtained by Alhassid et 
al. [1, Q . Thus the increase of the moment of inertia can 
be considered as another signature of the pairing phase 
transition together with the 5-shape of the heat capacity 
44 
42 



43 
41 



and the suppression of odd-even mass differences 
0. 

In the next step, we investigate the effect of the 
temperature-dependent moment of inertia on the spin 
distribution of nuclear levels for the case 56 Fe. For this 
purpose, we plot for the temperatures T = 0.5, 0.67, 1.0, 
and 2.0 MeV the spin distribution p(E, I) / p(E) from Eq. 

in Fig. O Moreover, we plot as function of temper- 
ature the first [i\ and second psi moments of this spin 
distribution where 



/'2 



(71) 
(72) 



and fj = p(E, I)/p(E) in Fig. [71 and finally, we plot again 
as function of temperature relative components of the 
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FIG. 7: (Color online) First (circles) and second (triangles) 
moments of the spin distribution as function of temperature 
for 56 Fe. Results from the SPA+RPA calculation (full sym- 
bols) are compared to a model using the rigid-body moment 
of inertia (open symbols) . 



spin distribution for the spins 7 = 0,2,4, and 6 h. Com- 
paring the SPA+RPA results to the spin-cutoff model 
based on a rigid-body moment of inertia, we observe from 
Figs. [6j B and [8] that good agreement between the two is 
achieved at high temperatures T > 2 MeV. Deviations 
are, however, significant for lower temperatures. In par- 
ticular, the 1 = component is much larger than all other 
/-components in the SPA+RPA calculation, whereas for 
the rigid-body spin-cutoff model, the largest component 
at low temperature is / = 2 while the 1 = component 
is quite small. Around T - 0.5 MeV in the SPA+RPA 
calculation, the 1 = component quickly decreases with 
increasing temperature whereas the 1 = 2 component 
increases drastically up to T = 0.7 MeV where it be- 
comes dominant, in agreement with the rigid-body spin- 
cutoff model. The 1 = 4 and 6 components also exhibit 
small magnitudes at low temperature and increase with 
increasing temperature. 

To conclude the discussion, we would like to make two 
remarks. Firstly, the spin distributions in Fig. [6] can be 
very well described by a phenomcnological model for the 
moment of inertia 



Cl 



1 + exp(c 2 T 2 ) 



C3, 



(73) 



with the parameters c\ = 15.0 h 2 /MeV, c 2 = 1.8 MeV~ 2 , 
and C3 = 1.5 fi, 2 /MeV This model of the moment of 
inertia includes the effect of the quenching of pairing 
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FIG. 8: (Color online) Spin components of level density as 
a function of temperature T for 56 Fe. Relative components 
for spins J = 0, 2, 4, and 6 from SPA+RPA calculations 
(upper panel) and from a spin-cutoff model using the rigid- 
body moment of inertia (lower panel). 



correlations. The three parameters of the model essen- 
tially determine the moment of inertia at zero tempera- 
ture (in the presence of pairing correlations), at very high 
temperature (in the absence of pairing correlations), and 
the critical temperature itself. Hence, the three parame- 
ters might be connected to collective E2 transition rates 
within rotational bands at low and high spin, respec- 
tively, and to odd-even mass differences via T c sw 0.57A. 
This would be interesting to investigate since a simple 
phenomcnological model of the moment of inertia based 
on experimental observables would be valuable for many 
applications of the nuclear level density such as Hauser- 
Feshbach-type cross-section calculations. Secondly, we 
note the presence of an S shape in the second moment of 
the spin distribution (see Fig.[7|). This is in analogy to the 
S shape of the heat capacity (see Fig. ^ which in itself 
is proportional to the second moment of the excitation- 
energy distribution for a given T. We would therefore 
propose this S shape of the second moment as another 
signature of the pairing phase transition, and it would be 
interesting to look for experimental evidence for it. 
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VI. CONCLUSION 

We have presented a treatment of nuclear rotation by 
the SPA+RPA method including thermal and quantal 
fluctuations. Using the collective coordinate method and 
applying a Marshalek-Weneser treatment, the RPA par- 
tition function is separated into two parts: an intrinsic 
RPA and a rotational part. We then can derive the spin- 
dependent level density p(E, I) in the saddle-point ap- 
proximation without the assumption of individual ran- 
dom spin vectors. The spin-cutoff parameter a can be 
identified with the moment of inertia J x of the gener- 
alized Belyaev formula. As an illustrative example, we 
have applied this method to 56 Fe and investigated the 
thermal properties including the moment of inertia with 
consideration of thermal and quantal fluctuations. We 
obtained the same conclusion as Alhassid et al. [3, Q, 
namely that the increase of the moment of inertia with 
increasing temperature is attributed to the suppression 
of pairing correlations. The increase of the moment of 
inertia can therefore be considered as one of the signa- 
tures of the pairing phase transition in parallel with the 



5-shape of the heat capacity and the suppression of odd- 
even mass differences. The spin distribution of the level 
density shows that the 1 = component is dominant at 
low temperatures due to the presence of pairing correla- 
tion; this component decreases with increasing tempera- 
ture. At high temperatures, the level density is governed 
by / ^ components. In this work, we assumed an 
axially symmetric nucleus with K = although our for- 
mulation can be easily extended to non-zero K quantum 
numbers. However, the K quantum number is thought 
to disappear due to if-mixing at some excitation energy 
above the yrast line. For this mixing, another degree of 
freedom, namely 7 deformation is important. It would 
be very interesting to investigate this issue in the future. 
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